Approximate vortex solution of Faddeev model 
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Through an Ansatz specifying the azimuthal-angle dependence of the 
solution, the static field equation for vortex of the Faddeev model is con- 
verted to an algebraic ordinary differential equation. An approximate 
analytic expression of the vortex solution is explored so that the energy 
per unit vortex length becomes as small as possible. It is observed that the 
minimum energy of vortex is approximately proportional to the integer 
which specifies the solution. 

PACS: 11.10.Lm,02.30,Ik,03.50-z 

1 Introduction 

Faddeev model PP was originally proposed as a model which might give rise to 
3-dimensional soliton solutions. Later it was discussed that the model might 
be an effective field theory describing the low energy behavior of the £{7(2) 
gauge field[l]. It can be regarded as a restricted version of the Skyrme model[3 
which may be an effective field theory for hadron dynamics. Because of the 
high nonlinearity of Faddeev and Skyrme models, the analytic structures of 
solutions of these models have not yet been clarified. The numerical solutions of 
these models, however, exhibit quite interesting soliton properties [3]- [pj: soliton 
solutions of the Skyrme model have tetrahedral structures while those of the 
Faddeev model have knot structures 0],[S]. 

The Faddeev model is a model concerning the real scalar fields 




(1) 
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satisfying 

3 

n 2 (x) = n{x) ■ n(x) = ^ n a (x)n a (x) = 1. (2) 

o=l 

It is defined by the Lagrangian density 

Cf(x) = c 2 l 2 (x) + gJ 4 (x), (3) 

h{x) = d^nix) ■ d»n{x), (4) 

i 4 (aO = -H„ v (x)H^(x), (5) 
H^ v {x) = n(x) ■ [d^nix) x <9„n(a;)] 

= e abc n a (x)d t ,n b (x)d u n c (x), (6) 

where c 2 and c 4 are constants. The static energy functional -Ef[«,] associated 
with Cf(x) is given by 

E F [n] = J dVe{x), (7) 
e(x) = c 2 e 2 (x) + c 4 e 4 (x), (8) 

3 3 

e2(x) = J2Y,[d^(x)r, (9) 

a— 1 i—1 

C4 (x)= ^[^(x)] 2 , (10) 

with a; = (xi, £2, 23) and c£V = dx\dx 2 dx%. By the stereo-graphic projection, 
the field n can be expressed by a complex function u as 

„_( u + u * -i{u-u*) \u\ 2 -l \ 

lH 2 + i' H 2 + i 'h 2 + J- ( j 

In terms of u, the energy densities e 2 and e 4 are given by 

£2 =(iTW (VM ' W) ' (12) 
64 - ~ 8 (1 + H 2 ) 4 • (13) 

The field equation can be rather simply expressed by q defined by 

q = XVu, (14) 

with 

u = Re^, R=\u\ (15) 

and 
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where 1\J C4/C2 of the dimension of length has been set equal to 1 and q is 
dimensionless. The static field equation can be written as [7] 



V -a + i/3 •£* = (), (17) 
where a. and (3 are a complex and a real 3- vectors given by 

a = q*-q* x(qxq*), (18) 
/3 = -(u*q - uq*) = BV$, B = (19) 

respectively. 

In this paper, we discuss some simple analytic solutions of the above equa- 
tion. We shall consider mainly the vortex solutions of the form R — R(p), $ = 
mcj) with (p, z, (j>) and m being the cylindrical coordinates and an integer, respec- 
tively. We adopt an approximate analytic solution containing two adjustable 
parameters. Its form is fixed by the compatibility with the boundary condition 
and the singularity-structure necessitated by the field equation. The parameters 
are fixed so that the energy per unit length of vortex, A m , becomes minimal. 
We find that the minimum of A m is proportional to m approximately. 
This paper is organized as follows. In Sec. 2, we discuss the m = case briefly. 
In Sec. 3, the above-mentioned approximate analytic solution for the m ^ case 
is discussed. Sec. 4 is devoted to summary. 

2 Some solutions in special cases 

We first consider the solutions of the form 

R = R(p,z), $ = $(0), (20) 
where (p, z, 4>) are cylindrical coordinates. Then we have 

q = Xc l * lR p e p + R z e z + ——eA , (21) 

a = e-' ; * (Ce p + De z + Fe ) , (22) 
B$' 

f3 = e , (23) 



D = 


[ x + 


F = 




—1 




P 


Y = 


2X 3 , 


= 


d<f>(<t>) 



p 

R 2 Y(<S>' 

7~ 



Rz, (25) 



[X + Y (R 2 p + Rl)] , (26) 

(27) 



d R(p,z) Rz = dR( P , z) 
dp dz 



u p — ~ ' Uz ~ aZ — ' \ z °) 
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where e p ,e z , and are orthonormal unit vectors satisfying e p x = e z , etc. 
We obtain = from Im{V ■a + i/3-a}=0. Then, taking the single- valuedness 
of u into account, we are led to 

$ = m<p, m : integer. (29) 

On the other hand, from Re{V • a + i/3 ■ ct}=0 and $ =m, we have 

V- (XGVR) + (G - l) R ~ 1 [l + 2X 2 (Vi?) 2 l = 0, (30) 
2R L 

2m 2 R 2 X 2 , s 

G = l + = . 31 

P 

In the m = case, it becomes Afarctani?] = 0, whose solution is given by 

°° / B \ 

arctani? = £ ( A„r" + J P„(cos0), (32) 



r = \J p 1 + z 2 , cos6> = -, (33) 
r 

where {A n , B n : n = 0, 1, 2, • • • } are constants and P n denotes the n-th Legendre 
polynomial. It should be noted that 

arctani? = ax lnp + a2 (ax, a-i : const.) (34) 

is a solution, too. These results indicate how Cho's solutionJS] u — co " st ' should 
be generalized. They are, however, trivial in the sense that the interference ef- 
fects between the ^(x) and £4(2;) in the Lagrangian disappear in these solutions. 
For m / 0, we see that the simpler cases such as R = fx{z),R = f%(r), R = 
fz{6), R = fn{rj), R = fz{C) are not allowed since G in Eq. (j3Tj) contains p, 
where (r,9,<fi) and (?7, are polar- and toroidal- coordinates, respectively. 
Only the vortex solution R = R(p), which is a solution of the 2-dimensional 
Faddeev model, is allowed. 



3 Approximate vortex solutions 

We hereafter consider the case R — R(p) and m 7^ 0. If we set 

R = tan-, a = — , (35) 

we obtain 

(a 2 + 2sin 2 44 + sin2£ ( + ( a - -sin 2 A ^ ~ 2m 2 sin2£ 
v ' da z \aa ) \ a J da 

Through a further change of variables 



0. (36) 



V(0 = -t^ = -[ w - 1 , C =^ (37) 



2R x 2 
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we obtain an algebraic differential equation 



cPV_ _ 1 / 1_ , J_ , 3 \ fdV^ 2 

1 1 \ dV 2m 2 V(V - 1) 



c-i c-vj d( ai-o 2 (v-c) 



= 0. (38) 



The solution of this equation can be explored in the following way. For a given 
Co which is different from and 1, we assume that the solution V(() near C = Co 
is given as 



nC) = E F «(C-Co) a+n , 0<Co<l, a<0. (39) 

n=0 

Then we are led to 

a = -2 (40) 

and 

Vr = (41) 
Co - 1 

v (4m 2 + Co)Vo-4Co(Co-l) 2 (Co-3) 

V2 = 12Co(Co - I) 2 ' (42) 

m 2 T/ + 2Cq 2 (Cq - I) 2 
V3 ~ 6 Co 2 (Co - I) 2 ' (43) 

and so on. It turns out that Vq and Co are arbitrary and Vi, V2, V3, ■ ■ ■ are fixed 
by them. 

Although V n can be determined order by order, it is difficult to conclude that 
the series (3.5) converges in some domain of C around Co- To obtain the radius 
of convergence of the series (3.5), if any, we must calculate V n for very large 
n, whose mathematical expression becomes quite complicated as n gets large. 
The differential equation (3.4) suggests that the behavior of V(C) at C = and 
1 (p = and 00) may be different from that at ordinary points < C < 1 
(0 < p < 00). It is also difficult to maintain the properties of V(Q near 
C = and 1 by an infinite series of the form (3.5). We therefore consider an 
approximate V(() which realizes the properties required by the lowest order 
analysis of (3.4) around C = Co 0j 1)> C = and C = 1- 

As was explained above, the behavior of V(C) in the neighborhood of Co is given 
by V(Q ~ const. (C — Co) -2 - We next consider the behavior of V(Q near C = 0. 
It turns out that the behavior V(C) ^const.C A is compatible with (3.4) only 
when A = 2m or A = —2m. We here recall that, with the assumption (2.10), 
only R = and R = 00 yield the vanishing energy density irrespective of the 
value of (f>. These values of R imply V(0) = 0, leading to A = 2|m|. In other 
words, if we require that the energy density is vanishing at p — 0, the differential 
equation (3.4) and the assumption (2.10) lead us to 

V(C) ~const.C 2|m| , C~0. (44) 
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Similarly, from the requirement that the energy density is vanishing at p = oo 
, we obtain 

^(C)~const.(C-l) 2H , C~l- (45) 

It is straightforward to see that the conditions V(0) = V(l) = correspond to 
the configurations n = (0,0,1) or n — (0,0,-1) at p — and oo. Therefore 
we are here considering the configurations interpolating these configurations. In 
the example considered below (Fig. 4), we obtain V(() connecting n — (0,0, 1) 
at p = and n = (0, 0, —1) at p — oo. 

We note that, in the analysis of the hedgehog Skyrmion, the solution of the 
differential equation 

d 2 W 1/1 1 t 3 \ fdW\ 2 

~d~q 2 ~ ~ 2 \W + rj-W + W-l) \~d~n) 
^1/1 1\ 1 



2 \ 77 — 1 77/ 7] — W 



was investigated[9] by a trial function 



dW _ W[(r) + 1)W - 2rj\ _ 
~d~^ ~ 2r] 2 (i]~ 1) 2 {W -if) ~ 

(46) 



W= - q ^- l) \ 9 >0, 1> P >0, (47) 

v-p 

which yielded a rather good value of the energy. 

In the following, we consider the approximate solutions V m (£) given by 
\r(l - r)] 2m 

V m (C) = -qm l \\ % , m>0, q m >0, l> Pm >0, (48) 

IS. Pm) 

which is compatible with the singularity-structure indicated by the differential 
equation and the boundary condition that the energy density vanishes at p — 
0, 00. 

The energy Ep corresponding to the assumption (2.10), which we denote by 
E m , is now given by 



E m = 87Tc 2 / A m dz, (49) 



1 



A m = I W m dC, (50) 

^[C(l-C)] 2m_1 







" 8 (c~p m ) 2 +q m m-o] 2m 

1 q m [C(l - C)] 2 " 1 ' 1 [C(l - + m(2( - IXC - p m f 

2 {(C-p m ) 2 + qnAC(l-C)} 2m } 2 

| 1 q 2 m C 4m - 2 (l C) 4m [C(l ~ + m(2C - IXC ~ Pm)? _ (51) 
8 {(C-p m ) 2 + q m [C(l-C)} 2m } 3 
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We fix p m and q m so that they minimize A m . For m = 1 ~ 4, they are given as 



m = 1 
m = 2 
m = 3 
m = 4 



pi = 0.83, qi = 4.0, Ai = 1.14 

p 2 = 0.61, g 2 = 2.3 x 10, A 2 = 2.23 

p 3 = 0.56, q 3 = 1.57 x 10 2 , A 3 = 3.28 

p 4 = 0.52, q 4 = 1.36 x 10 3 , A 4 = 4.34. 



It seems that p m approaches to 0.5 when m becomes large. For the con- 
venience of numerical estimation, we fix p m as 0.5 for m larger than 4 and 
determine q m so as to minimize A m . The result is given as follows: 





m — 5 




g 5 = 1.35 x 10 4 


A b = 5.42 




m — 6 




q 6 = 1.45 x 10 5 


A 6 = 6.49 




m = 7 




q 7 = 1.66 x 10 6 


A 7 = 7.57 


< 


m = 8 




q s = 1.96 x 10 7 


A s = 8.65 




m = 9 




q 9 = 2.41 x 10 8 


A 9 = 9.74 




to = 10 : 


gio = 3.06 x 10 9 


A w = 10.8. 


The 


results for ^4 


„ and qm+i/q-m are given in Fig.l and Fig. 2, respectively. We 



are then led to the approximate formulae 



A m — 



bm, a = 0.052, b = 1.076 



and 



c (1 - d e- fm ) 



14.0, d= 0.769, / = 0.234. 



(52) 



(53) 



Since the infinite product Ilm=o (l — °^ c _ -^ m ) = B converges to 0.00592, we 
have log q m ~ m log c + log B for large m. From the fact that a ~ and 6 ~ Ai, 
it looks like that there exists an atom-like object Q with energy b and unit m 
and that the assembly of m Q's constitutes the configuration with to. 

With the original variables, the energy density K (p) defined by A\ — J °° K{p)dp 
and R(p) for m = 1 are given in Fig. 3 and Fig. 4, respectively. 




Figure 1: A m vs. to 
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4 Summary 

We have examined the vortex solutions of the Faddeev model. Through the 
change of variables, the field equation was converted to an algebraic differen- 
tial equation containing an integer-parameter to. Its approximate solution was 
parameterized with the aid of two parameters p m and q m . They were fixed so 
that the energy per unit length of vortex, A m , became minimal. It was observed 
that the minimum of A m is proportional to to approximately. The numerical 
analysis of the Faddeev model made so far clarified the knot structure of the 
genuine three-dimensional solitons. We hope that our two-dimensional analysis 
might be helpful for the understanding of the latter since knot-solitons may be 
regarded as the bended and twisted vortices. 
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